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Abstract —In this paper, a statistical model for the coupling 
of electromagnetic radiation into enclosures through apertures is 
presented. The model gives a unified picture bridging determinis¬ 
tic theories of aperture radiation, and statistical models necessary 
for capturing the properties of irregular shaped enclosures. A 
Monte Carlo technique based on random matrix theory is used 
to predict and study the power transmitted through the aperture 
into the enclosure. Universal behavior of the net power entering 
the aperture is found. Results are of interest for predicting the 
coupling of external radiation through openings in irregular 
enclosures and reverberation chambers. 

Index Terms —Statistical Electromagnetics, Cavities, Aperture 
Coupling, Admittance Matrix, Chaos, Reverberation Chamber. 


I. Introduction 

T he coupling of electromagnetic radiation into enclosures 
or cavities through apertures both electrically small 12 
and large 1^ has attracted the interest of electromagnetic 
community for many years Ei-ii. Full solutions of this prob¬ 
lem are particularly complicated because of the mathematical 
complexity in the solution of the boundary-value problem 
and because of the sensitivity of the solution to the detail 
of the enclosure’s dimensions, content, and the frequency 
spectrum of the excitation. These difficulties have motivated 
the formulation of a statistical description (known as the 
random coupling model, RCM IS), 16], Q) of the excitation 
of cavities. The model predicts the properties of the linear 
relation between voltages and currents at ports in the cavity, 
when the ports are treated as electrically small antennas. 

In this paper, we formulate and investigate the RCM as 
it applies to cases in which the ports are apertures in cavity 
walls. The aperture is assumed to be illuminated on one side 
by a plane electromagnetic wave. We then distinguish between 
the radiation problem, where the aperture radiates into free 
space, and the cavity problem, where the aperture radiates into 
a closed electromagnetic (EM) environment. The solution of 
the problem in the cavity case is then given in terms of the 
free-space solution and a fluctuation matrix based on random 
matrix theory (RMT). Thus, there is a clear separation between 
the system specihc aspects of the aperture, in terms of the 
radiation admittance, and the cavity in terms of the fluctuation 
matrix. 
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We illustrate our method by focusing on the problem of an 
electrically narrow aperture, for which the radiation admittance 
can be easily calculated numerically. We then generalize our 
result to the interesting case where a resonant mode of the 
aperture is excited. In this case the statistical properties of 
the aperture-cavity system can be given in a general universal 
form. 

Our results build on previous work on apertures. In partic¬ 
ular, rectangular apertures have a very long research tradition 
in EM theory 0, and continue to be a topic of interest 
0. The first self-consistent treatments have been carried out 
by Bethe IT], Bouwkamp 0, and, later, by Schwinger in 
aperture scattering Eo], and Roberts M- Subsequent work 
on apertures is due to Ishimaru 112, Cockrell E], Harrington 
0, 0, M, 113, and Ramat-Sahmii E], 13, among other 
investigators. 

Our results are of interest for the physical characterization of 
the radiation coupled into complex cavities such as reverbera¬ 
tion chambers (RC) - which is known to be an extremely com¬ 
plicated problem even challenging classical EMC techniques 
113, M, M - for understanding interference in metallic 
enclosures, as well as for modeling and predicting radiated 
emissions in complicated environments. 

The paper is organized as follows. In Sec. II. we introduce 
the general model for the cavity backed aperture, and we 
describe the way the RCM models the cavity. In Sec. III. we 
apply the formulation of Sec. II. to large aspect ratio, rectangu¬ 
lar apertures; evaluating the elements of the admittance matrix 
and computing the power entering a cavity with a rectangular 
aperture. In this section we also develop a simple formula for 
the power entering a low loss cavity with isolated resonances. 
Sec. IV. describes an extension of the model that accounts for 
the coupling of power through an aperture, into a cavity, and 
to an antenna in the cavity. Simple formulas for the high-loss 
case, low-loss, isolated resonance are developed. 

H. Random coupling model for apertures 

The random coupling model was originally formulated to 
model the impedance matrix of quasi-two-planar cavities with 
single 0, and multiple l6] point-like ports 0. In this section, 
we develop the model for three-dimensional irregular enclo¬ 
sures excited through apertures ll20l . As depicted in Fig. [T] 
we consider a cavity with a planar aperture in its wall through 
which our complex electromagnetic system is accessed from 
the outside. The size and shape of the aperture are important 
in our model, their specification constitutes “system specific” 
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Fig. 1. Geometry of a complex 3D cavity with volume boundary 

gycav^ exhibiting wave-chaotic ray trajectories in the semiclassical limit: 
energy enters the cavity through an electrically large aperture dV^P. Also 
shown is a point-like port to be considered in section IIVI 


information that is needed to implement the model. The cavity 
that we consider in our studies is an electrically large enclosure 
with an irregular geometry. The irregularity is assumed to be 
such that ray trajectories within the cavity are chaotic through¬ 
out Il2n . Generally speaking, typical cavities have this feature; 
particularly those with curved walls and/or with contents that 
scatter radiation in multiple directions. The consequences of 
assuming trajectories are chaotic is that the spectra of modes of 
the cavity have universal statistical properties that are modeled 
by random matrix theory (RMT) Il22ll . Experiments have been 
carried out to test the predictions of RCM for quasi-two- 
dimensional complex enclosures coupled through electrically 
small (point-like) ports ||23, ll24l . ||25l. Experiments have also 
been carried out in three-dimensional enclosures, and with 
electrically-large ports ||26l, IIZTl . 

When apertures are considered, the port approximation 
invoked in the original derivation of the RCM is no longer 
valid, hence we need to consider the field distribution on the 
aperture surface. We consider a planar aperture, i.e., one that 
is not subject to boundary curvature. This is consistent with 
real-world situations such as 3D reverberation chambers (RCs) 
128], where the aperture is generally in a planar boundary. 

We first treat the aperture as if it existed in a metal plate 
separating two infinite half spaces. We refer to this as the “free- 
space” situation or “radiation” case. Suppose the port is treated 
as an aperture in a planar conductor whose surface normal n is 
parallel to the z-axis. The components of the fields transverse 
to z in the aperture can be expressed as a superposition of 
modes (an example is the set of modes of a waveguide with 
the same cross sectional shape as the aperture) 


Et = Y, Vs^g (Xi) , 

S 

(1) 

Ht =^Ignx tg (xj^) , 

(2) 


S 


where Cs is the basis mode, (having only transverse fields) 
normalized such that J d^x± |es|^ = and n is the normal, 
which we take to be in the z-direction. In the radiation case we 
solve MaxwelTs equations in the half space z > 0 subject to 
the boundary conditions that Et = 0 on the conducting plane 
except at the aperture where it is given by O- We do this 


by removing the conducting plane, adding a magnetic surface 
current density 25{z)n x Et to Earaday's law, and solving 
MaxwelTs equations in the whole space — oo < z < oo in the 
Eourier domain. Eor this problem, the transverse components 
of the electric field are odd functions of z with a jump equal 
to twice O at the location of the aperture; thus satisfying 
the boundary condition for the half space problem. We then 
evaluate the transverse components of the magnetic field on 
the plane z = 0, and project them on to the basis fix Cg (xj^) 
at the aperture to find the magnetic field amplitudes Ig in 
(|2]l. The result is a matrix relation between the magnetic field 
amplitudes and the electric field amplitudes in the aperture, 

= , (3) 


where fco = w/c, uj is the frequency of excitation and we 
have adopted the phasor convention exp {—iujt). Here, the 
radiation admittance matrix is determined from the 

SS 

Eourier transform solution for the fields, and is given in terms 
of a three dimensional integral over wave numbers. 
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where the dyadic tensor 


Ayrad 


kl 


( k"^ - k\ {kl - k"^) k\ \ 

( j 

(k X n) (k X n) 

fcl ’ 


(5) 


is responsible for coupling two arbitrary modes of the aperture, 
and 

Cs = / d^xj. exp(-ikj. -x^jcs , (6) 

Japerture 


is the Eourier transform of the aperture mode. The elements of 
the radiation admittance are complex quantities. The residue 
at the pole fc = /cg in (|4|i gives the radiation conductance 



where is the two dimensional solid angle of the wave vector 
k to be integrated over dr, and where there appears a modified 
dyadic tensor 


^Qro.d = [(k_Lk_L) /fcy + [(k X n) (k X n) /k]^ . (8) 


The radiation conductance is frequency dependent through /cq. 
We note that there is an implicit fcg dependence through the 
Eourier transforms of the aperture modes, and we set |k| = feg 
in ® and The remaining part of (|4]i, gives the radiation 
susceptance. Part of this can be expressed as a principle part 
integral of the radiation conductance. However, there is an 
additional inductive contribution (Y oc kg^) (which we term 
the magnetostatic conductance) coming from the last term 
in the parentheses in (|5]) that contains a factor that cancels 
the resonant denominator in (|4]l. The reactive response of the 
aperture can be expressed in terms of the Cauchy principal 
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value of the radiation conductance ( 01 , and the previously 
mentioned inductive contribution, yielding 

B:f{ko)=pJ^ G7/ik) + B2f, (9) 

where stands for magnetostatic conductance, defined as 

and where Ab^.^ = [(k x fi) (k x n)] /k'^. 

We now repeat the process, but assume that the aperture 
opens into a cavity rather than an infinite half space, the 
radiation admittance ( 0 i will be replaced by a cavity admit¬ 
tance. In appendix it is shown that under the assumptions 
that the eigenmodes of the closed cavity can be replaced 
by superpositions of random plane waves (Berry’s hypothesis 
ISU, Gol), and the spectrum of the cavity eigenmodes can 
be replaced by one corresponding to a random matrix from 
the gaussian orthogonal ensemble (GOE) ll22ll . the statistical 
properties of the cavity admittance can be represented as 


ycav ^ 


rad 


[g 


^adl 




'adl 


( 11 ) 


and the matrix is a universal fluctuation matrix defined as 


= TT /Co — /Cn -f ia 

n 


( 12 ) 


where ICn is a set of eigenvalues of a random matrix drawn 
from the GOE Ii22^ . These matrices have positive and negative 
eigenvalues that fall in a symmetric range about zero. The 
mean spacing of the eigenvalues near zero can be adjusted 
by scaling the size of the matrix elements. Here we assume 
that this has been done such that the mean spacing of the ICn 
near zero is unity. The quantity /Co represents the deviation of 
the excitation frequency ojq from some reference value ujref 
placed in the band of interest 


_ ^0 ^ref 

Alu 


(13) 


where Aw is the mean spacing between the resonant fre¬ 
quencies of modes of the actual cavity in the range of uJref- 
The resonant frequencies of a particular realization are thus 
given by w„ = oJref + Aw/C„. The vector is composed 
of independent, zero mean, unit variance Gaussian random 
variables. The quantity a describes the average loss in the 
cavity and is related to the finesse parameter IF, widely used 
in optical cavities and photonic lattices ED 


x—1 

a = IF = ——— . 

2QAw 


(14) 


Thus, according to the RCM, two quantities are required to 
specify the properties of the cavity: the Q width ujq/Q and 
the mean spacing Aw between nearest neighbor resonances 
(cavity modes). 

The quantity a determines whether the cavity is in the high 
loss {a > 1) or low loss (a < 1) regime. More precisely, 
it establishes how much the individual resonances overlap. 
Basically, if a > 1 the Q-width of a resonance, ujqIQ, exceeds 
the spacing between resonances Aw. A feature of chaotic 


cavities is that each mode in a given frequency band has 
essentially the same (J-width. This is a consequence of the 
ergodic nature of the ray trajectories that underlie the mode 
structure. Each mode has wave energy distributed throughout 
the cavity such that losses are approximately the same for each 
mode. 

We now describe how relations (0i and (ED are used to 
determine the coupling of our system to external radiation. 
Specifically, we consider the framework of Eig. 0 where the 
aperture is illuminated by a plane wave incident with a wave 
vector k“° and polarization of magnetic field that is 
perpendicular to We again consider the aperture to be 
an opening in a conducting plane at z = 0 , with radiation 
incident from z = —oo. We imagine writing the fields for 
z < 0 as the sum of the incident wave, the wave that would 
be specularly reflected from an infinite planar surface, and 
a set of outgoing waves associated with the presence of the 
aperture. The incident and specularly reflected waves combine 
to produce zero tangential electric field on the plane z = 0 . 
Thus, the outgoing waves for z < 0 associated with the 
aperture can be expressed in terms of the electric fields in 
the aperture just as the outgoing waves for z > 0 can, the 
two cases being mirror images. So, relation E) continues to 
represent the tangential electric fields in the plane z = 0. Eor 
the magnetic field we have separate expansions for z > 0 
and z < 0. The electric field amplitudes are then determined 
by the condition that the magnetic fields are continuous in 
the aperture at z = 0. Eor z = 0-(- we have Hf = 
If h X is (x_l), with P" = ■ V, where Yf" is either 

the radiation admittance matrix ■01 or the cavity admittance 
matrix ETTi depending on the circumstance. Eor z = 0 “ we 
have Ifh x Cs ) -f 2 h*"'^ exp [/k™° • x_l] , where 

P = —Yf°‘‘^ ■ V (the minus sign accounts for the mirror 
symmetry) and the factor of two multiplying the incident 
field comes from the addition of the incident and specularly 
reflected magnetic fields. Projecting the two magnetic held 
expressions on the aperture basis, and equating the amplitudes 
gives 

+Yf‘'‘^) ■Y = 2P\ (15) 

where 

7 ™'= = -n • e, (-k™^) x (16) 

and Cs is the Eourier transform of the aperture electric held 
and is defined in (01. Equation ESl) can be inverted to And the 
vector of voltages 

y = + • 2/“% (17) 

and then the net power passing through the aperture is given 
by 

f, = r) , ( 18 ) 

where Y^ can be either a radiation (free-space) admittance 
@ or a cavity admittance matrix ED- 

III. Rectangular apertures 

We consider rectangular apertures and select a basis for 
representation of the tangential fields in the aperture in (01 and 
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Fig. 2. Geometry of a 2D regular (narrow) aperture illuminated by an external 
plane-wave, radiating in free-space. 


0. One choice for the basis are the modes of a waveguide 
with a rectangular cross section. These can be written either 
as a sum of TE and TM modes, or simply as a Fourier 
representation of the individual Cartesian field components 

ED Ga. 

In EMC studies, narrow apertures are often considered, as 
they frequently occur in practical EM scenarios. The aperture 
of Fig.|2]is elongated and thin: it has only one electrically large 
dimension, and the field component e(„ o.x) ~ 0- Hence, in the 
particular case of a rectangular aperture with L ~ A ^ VE and 
W ^ 0, the field will be dominated by TEno modes 

e. =» ™|fc,(^+V2)| . _ 

for \y\ < Wl2, where N = ^JYJLW. Once the aperture 
field basis is specified, the procedure for calculating the cavity 
admittance matrix (|4|i is given by the following steps ll34l . 
First, calculate the Fourier transform Cg (01. Second, use Cg to 
calculate the radiation conductance 0. Third, use Cg to cal¬ 
culate the magnetostatic susceptance (fTol i. Fourth, use above 
quantities to form the radiation susceptance 

ss ss 

Finally, use the so-formed radiation admittance to generate the 
cavity admittance (HB, similar to 0, m- The conductance 
and susceptance have the property that off diagonal terms 
vanish if n is odd and n is even, and vice versa, due to 
the even and odd parity in x of the basis modes. Further, in 
the limit of high frequency, ko ^ kn, the components 
have the limits 


G 


rad 

nn 



n = n , 
n ^ n , 


( 20 ) 


Thus, at high frequencies the admittance matrix is 

diagonal and equal to the free space conductance as would 
be expected when the radiation wavelength is much smaller 
than the aperture size. 

We have evaluated the elements of the radiation conductance 
matrix as functions of frequency for a rectangular aperture 
with dimensions L = 25 cm xW = 2 cm. Plots of these 



Fig. 3. Diagonal radiation conductance for n = 1... 5, for a narrow aperture 
L = 25 cm X W = 2 cm. In the inset: low frequency (cutoff) detail of the 
diagonal elements with an off-diagonal element. 



Frequency [GHz] 

Fig. 4. Close-up of the diagonal radiation conductance in Fig. [3] low 
frequency (cutoff) detail of the diagonal elements with an off-diagonal element 
(dashed line). 


appear in Figs.0and0 For the diagonal elements the conduc¬ 
tance increases nonmonotonically from zero, and asymptotes 
to its free-space value. For the off-diagonal components (not 
shown) the conductance first rises and then falls to zero with 
increasing frequency. Combining numerical evaluations of the 
two contributions in 0 yields the total susceptance. This is 
plotted for several diagonal elements in Fig. 0 which have 
the feature that for some elements the susceptance passes 
through zero at a particular frequency. At these frequencies the 
conductances (Fig. |4| also tend to be small. This indicates the 
presence of resonant modes of the aperture. At the frequencies 
of these modes the aperture will allow more power to pass 
through the plane of the conductor than would expected based 
on the aperture area. Our next step is to evaluate the power 
passing through the aperture when it is illuminated by a plane 
wave. This involves solving (fTTl l for the vector of voltages V, 
and inserting these in (fTSl l for the power transmitted through 
the aperture. This will be done a number of times; first with 
= YJ°‘^ to determine the power passing through the 
aperture in the radiation case; then again with Y> = 
to determine the power through the aperture when it is backed 
by a cavity. In the cavity case a number of realizations of the 
cavity admittance matrix will be considered, modeling cavities 
with different distributions of resonant modes. 

Figure 0 shows the frequency behavior of the transmitted 
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Fig. 5. Diagonal radiation susceptance for n = 1... 5, for a narrow aperture 
L = 25 cm X W = 2 cm. 



Fig. 6. Frequency behavior of the net power transmitted by an aperture 
A = 0.25 m x0.02 m in free space, = 1,..., 7 modes, for an external 
plane wave of = 1 mA/m, at oblique incidence <j) = 0,0 = 7r/4, and 
polarization = 7r/2. 


power defined in ( fTSl l. with for an oblique plane 

wave incident with angles 4>p = n/2, (j) = 0, and 9 — 7r/4, 
according to the coordinate frame of Fig. ID In Fig. |6] the 
net power is parameterized by the number of aperture modes 
included in the calculation, ranging from = 1 to = 7. 
As expected, the higher the frequency the greater the number 
of modes required to achieve an accurate prediction of the 
transmitted power. Interestingly, we notice the presence of a 
sharp peak at 600 MHz, which is the slit resonance frequency, 
i.e., Ja = c/(2L), and other broader resonances located at 
n Ja- At normal incidence, the peak of the resonance next to 
the sharp peak is reduced of a factor of 10. This is confirmed 
in previous studies based on the transmission line model of a 
narrow aperture 1^ . Having computed the radiation conduc¬ 
tance and susceptance for a narrow slit of dimensions 25 cm 
X 2 cm, we can investigate the effect of a wave-chaotic cavity 
backing the aperture by replacing by and exploiting 
the statistical model, We first calculate distributions of 
the cavity admittance elements from the RCM by using a 
Monte Carlo technique and the radiation admittance of a 
rectangular aperture. Numerical calculations of dUi and Monte 
Carlo simulation of (Ell allow for generating an ensemble 
of cavity admittances of the form (fTTT i. In particular, we use 
the statistical method described in Refs. 0, and 0, with 
N = 7 aperture modes, M = 600 cavity modes, and a loss 
factor of a = 6, simulating a chaotic cavity with high losses. 



Fig. 7. Comparison between the radiation (red dashed line), and the average 
(red squares) net power transmitted by an aperture of dimensions A = 0.25 
m x0.02 m, in the frequency range from 1 GHz to 1.2 GHz for an external 
incidence of = 1 A/m, direction of incidence (p = 0,0 = 7r/4, 

polarization cpp = 7r/2. Thin solid (black), dashed (blue), and dash-dotted 
(puiple) lines are reported to show 3 independent cavity realizations as 
generated through a Monte Carlo method. The thick solid (red) line indicates 
the ensemble average of the port power over 800 cavity realizations. The thick 
dashed (red) line indicates the power radiated from the aperture in free-space. 
Each cavity response is given by the superposition of 600 ergodic eigenmodes. 
The chaotic cavity is modeled with a = 1.0. 


to create the bare fluctuation matrix (El- Here, we repeat 
the simulation of (El 800 times to create an ensemble of 
fluctuation matrices. By virtue of its construction, the average 
cavity admittance equals the radiation matrix = yrad 

Further, in the high loss limit (a ^ 1) fluctuations in the 
cavity matrix become small and For finite 

losses the character of the fluctuations in the elements of the 
cavity admittance matrix change from Lorenzian at low loss 
{a <C 1) to Gaussian at high loss (a ^ 1). 

We now consider the net power coupled through an aperture 
that is backed by a wave-chaotic cavity. This involves evalu¬ 
ating (flTl i and ( fTSl l with and with evaluated 

according to (El- When this is done the frequency dependence 
of the net power acquires structure that is dependent on the 
density of modes in the cavity. This is illustrated in Figs. 
|7] and [8] where net power through the 0.25 m by 0.02 m 
rectangular aperture is plotted versus frequency in the range 
1.0 < / < 1.2 GHz for a few realizations of the fluctuating 
admittance matrix. In both cases the reference frequency in 
(El is set at fref = i^ref I{“^Tr) = 1.1 GHz, and the mean 
spacing between modes is set to be A/ = Aa;/( 27 r) = 10 
MHz. Figure |7] corresponds to a moderate loss case {a = 1.0), 
and Fig. 0 to a low loss case {a = 0.1). Clearly, in the 
low loss case the peaks in power associated with different 
resonances are more distinct and extend to higher power. Also 
plotted in Figs. |7]and|8]are the average over 800 realizations 
of the power coupled into the cavity and the power coupled 
through the aperture in the radiation case. Notice that both of 
these curves are smooth functions of frequency and that in the 
moderate loss case of Fig. |7] the power averaged over many 
realizations (Pcav) is only slightly less than transmitted power 
in the radiation case. In the low loss case of Fig.[8]the average 
power coupled into the cavity is about half the value of the 
radiation case. 

The features of the coupled power in Fig. |7] and 0 can be 
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Fig. 8. Comparison between the radiation (red dashed line), and the average 
(red squares) net power transmitted by an aperture of dimensions A = 0.25 
m x0.02 m, in the frequency range from 1 GHz to 1.2 GHz for an external 
incidence of = 1 A/m, direction of incidence cp = 0,0 = 7r/4, 

polarization (pp = 7r/2. Thin solid (black), dashed (blue), and dash-dotted 
(purple) lines are reported to show 3 independent cavity realizations as 
generated through a Monte Carlo method. The thick solid (red) line indicates 
the ensemble average of the port power over 800 cavity realizations. The thick 
dashed (red) line indicates the power radiated from the aperture in free-space. 
Each cavity response is given by the superposition of 600 ergodic eigenmodes. 
The chaotic cavity is modeled with ol = 0.1. 


captured by a simple model. We consider the frequency range 
in Fig. |6] where the resonances in the aperture are isolated, 
roughly speaking / < 3 GHz. In this case, in a given narrow 
frequency range the aperture fields are dominated by a single 
mode, and we can replace the matrix equations (I15II18I) with 
their scalar version 


P^=2 SR 


_I!_\ 

|yrad_^y>|2j ’ 


( 21 ) 


where or depending on whether a cavity is 

present, and y“" = 

with 

^ = - E F- r ■ (22) 

TT ^ ALo — l^n’ + *0^ 

n 

Here /Co, /C„', and a are defined as before following (fl^ and 
the w^i are a set of independent and identically distributed 
Gaussian random variables with zero mean and unit variance. 
The quantities and are properties of the aperture, 
and will be discussed below. 

The power transmitted through the aperture in the absence 
of a cavity is given by 


where 


A. = 


2Qa (w - loa) 


UJA 


(25) 


and the effective quality factor for the aperture is given by 


n 

2Gr-^UA) du ^ ^ 

When a cavity backs the aperture, becomes a 

random frequency dependent function through the variable 
^ (w). This quantity is frequency dependent through the de¬ 
nominators in (l22l) and is random due to the random vectors of 
coupling coefficients Wn and eigenvalues /C„. The frequency 
scale for variation of ^ (w) is determined by the frequency 
spacing of modes of the cavity. In the typical case the 
frequency spacing of cavity modes is much smaller than that 
of aperture modes, as depicted in Figs. [T] and The behavior 
of the coupled power as a function of frequency will thus 
follow the envelope of the radiation case, with fluctuations on 
the frequency scale of the separation between cavity modes. 

This behavior can be captured in the simple mode if we 
assume the frequency is close to one of the poles (n = n) of 
(l22l) . Specifically we write 


£,=ib + 


TT (/Co — /C„ + ia) ’ 


(27) 


where the term i 6 is in the form of (l22l l with the n = n term 
removed and /Cq = /C„. Since a is assumed to be small it can 
be neglected in b (making b purely real) whereas it is retained 
in the second term in (l27l) since we consider frequencies such 
that /Co — ICn is small and comparable to a. The statistical 
properties of the sum given by b were detailed by Hart et al. 
IIJTII . If we express b in the form b = tanijj then the PDF 
of ip is cos'0/2 IIJTI . For now, since we are focusing on the 
behavior of individual realizations, we leave the value of b 
unspecified. 

Having assumed the cavity response is dominated by a 
single resonance we can manipulate (l2Tl l into a general form 


Pt{uj) = PA- 


AaaA 


|/Co ~ /C„ + / (a + a a) \ 
where K^ = Kn — 2A ua, A = Aa + b/2, and 


(28) 


a A = 


TT |1 + 2/A' 


2 ■ 


(29) 


Pt = 


jinc 

“T" 


^rad 


\Qrad iprad^ 


2 • 


(23) 


The aperture resonance occurs at a frequency oja where 
j^rad (^^) = 0, corresponding to a peak in Fig.|6] and the cou¬ 
pled power at the peak is Pa — Pt (oja) = | /(2G’'“'^), 

with G'’“‘^ = G^°'^{uja)- We can then express the frequency 
dependence of the power coupled through the aperture for 
frequencies near uja in the form of a Lorenzian resonance 
function 


Ptiu^) 


Pa 

1 + 


(24) 


Here, the variables have the following interpretations; ua is 
the “external” loss factor describing the damping of the cavity 
mode due to the aperture. Note, it is added to the internal loss 
factor a in the denominator of (l28T l. It is a statistical quantity, 
mainly through the Gaussian random variables Wn- This is 
responsible for variation in the height of the peaks in Fig. [8] 
The quantity A = A^i -1-6/2 represents the modification of 
the aperture resonance function (l25l) by the reactive fields of 
the nonresonant cavity modes (5/2). Note that it affects the 
external damping factor ua, which is largest when A = 0, 
i.e., when oj is near the aperture resonant frequency. Finally, 
/C„ determines the shifted cavity mode frequency. That is. 
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considered the case of fields excited by a current distribution 
of the form 1^ 

= (31) 

p 

where Up (x) is a set of basis functions used to represent the 
current distribution in terms of a set of amplitudes Ip which we 
call the port currents. The corresponding port voltages were 
dehned in Ref. as 

Vp = - J (fxup{x)-Eix) , (32) 


Fig. 9. Universal aperture factor on A a, parameterized by a. 

using definition (fT3T i the resonant cavity mode frequency, /Co = 
/C„, becomes 

uj = ujn — 2 AujA aA ■ (30) 

Equation (l28l l implies that the power coupled into the cavity 
at frequency w is bounded above by the power that can be 
transmitted through the aperture at the aperture resonance ioa- 
These powers are equal Pt (w) = Pa if the mode is resonant, 
/Co = /C„, and the cavity is critically coupled, a = aA- 
Note, however, that the coupled power in the cavity case (l28l l 
can exceed the radiation case (l24l i at the same frequencies if 
the aperture is off resonance A (w) ^ 0. This is evident at 
the peaks of the coupled power in Fig. Basically, what is 
happening is the cavity susceptance, which alternates in sign 
as frequency varies on the scale of the cavity modes, cancels 
the aperture susceptance, thus making the aperture resonant 
for frequencies away from the natural resonance uja- 

When the cavity loss parameter is small as in Fig. |8] or 
(l28l l there are large variations in coupled power as frequency 
is varied. A broad band signal would average over these 
variations. We can treat this case by computing the power 
coupled through the aperture averaged over realizations of the 
random variable ^ defined in (l22T l. Such averages are shown 
in Fig. Q and [8] based on Monte Carlo evaluations of the full 
system (HI. We perform this average in our simple model. 
A plot of a numerical evaluation of (Pt) /Pa from (l28l l as a 
function of A for different loss parameters a appears in Fig. 
|9] Interestingly, a loss parameter that is as large as a = 1 is 
sufficient to make the average power entering a cavity 90% of 
that passing through an unbacked aperture. 


and as a result the power entering the cavity through the ports 

is pp = ^{j:^Vp* Ip}/2. 

In analogy to our treatment of the aperture, we consider 
two cases; one in which the current distribution radiates 
into free-space and one in which the curi'ent distribution 
radiates into a cavity. The linear relationship between the port 
voltages Vp and the port currents Ip is then characterized by 
impedance matrices and depending on the case 

X pp pp X 

under consideration, wnere 

= (33) 

p 

For the radiation case it was shown ||38]| 



(Pk iko 
{2TTf kl - F 


Up-A^’V 


• (34) 


where Up {k) is the Fourier transform of the basis function 
Up (x), fco = ^^/c, and the dyadic Az is given by 


Az = 


\kr — kk 

P 


kk 

U2l2 

hj fVQ 


{kl - p 


(35) 


The radiation impedance matrix can be decomposed Z'', = 

^ ^ nn 


PP 


R, + iX^, where R^, is the residue from the pole at 

pp pp pp 

k = Iq in ( 13^ (the radiation resistance) 


(ko) = SR f Z^^f) 

pp k ^\ pp J 


kldflk 




lfc2 - kk 

P 


and is the reactive contribution. 

pp 

For the cavity case it was shown 




radi 


^ \M 


.radi 


• U / 

P 


(36) 


(37) 


IV. Cavity with both ports and apertures 

In the previous section we determined the properties of the 
power entering a cavity through an aperture. The contents of 
the cavity were treated as a distributed loss characterized by 
a single parameter a. We will now extend this model to treat 
the case in which we identify a second port to the cavity that 
we will treat as an electrically small antenna. This port could 
be an actual port where a connection is made to the outside 
world as illustrated in Fig. [T] or the port could represent the 
pin of a circuit element on which the voltage is of interest. 

We consider the conhguration of Fig. [T] where we have 
the joint presence of apertures and ports. We have previously 


where the fluctuating matrix ^ is defined in (fT^ . 

We are now in a position to describe a statistical model for 
a cavity including both an aperture and a localized current 
distribution. In this case we construct an input column vector 
(j) that consists of the aperture voltages and port currents and 
an output vector ijj that consists of the aperture currents and 
port voltages 




Ea 

Ip 


(38) 


and 


^ = 


Ia 

Yp 


(39) 

















where Va.p are the aperture and port voltages, and Ia^p are 
the aperture and port currents. These are then related by a 
hybrid matrix T, 4’ =T.‘ where 

X = t Im (g) + [Y] -i - [Z] ■ (40) 

Here the matrices U and V_ are block diagonal, viz.. 


U = 


^yrad 

“o 



(41) 


and V = ^ [^]. The dimension of U and V is {Ns + Np) x 
{Ns + Np), where Np is the number of port currents and Ns 
is the number of aperture voltages. Here we have assumed that 
the ports and apertures are sufficiently separated such that the 
off diagonal terms in U, describing the direct excitation of 
port voltages by aperture voltages, and aperture currents by 
port currents, are approximately zero. This assumption can be 
released in case of direct illumination between, or proximity 
of, aperture and port through the short-orbit correction of the 
RCM, which in this setting needs to be extended to cope with 
vector electromagnetic fields llJTll . 1391. In the simples case 
we can take the square root of V 


[^] 


1/2 


jradj1/2 
0 


0 

rad] 1/2 


(42) 


At this point we specialize consideration to the case of a port 
that is an electrically small antenna characterized by a single 
current Ip and single basis function in OTl l. The matrices T, 
U, and V then have dimension {Na + 1) x {Na + 1) and the 
matrix relations for the voltages and currents are more clearly 
expressed when separated into Na aperture equations and one 
port equation. For the aperture equations we find that (flST l is 
replaced by 

(r™"+ ■Ya+-Up ■ Ip = 2/“° ■ 

(43) 

For the port we assume the small antenna drives a load with 
impedance Zp such that Vp = —Zpip. Then the port equation 
becomes 


{Zp + iAP ■ -Zt = 0 , (44) 


where ^ap is a column vector defined similarly to the Na x 
Na matrix ^ in (fT^ . Finally, the scalar cavity impedance is 
defined in accord with dJTl i 

Zcav ^ -j^rad j^rad^ ^ ( 45 ^ 


with ^ being a scalar version of the fluctuation matrix. If we 
assume the aperture voltages are known, dSll can be solved 
for the current in the port, 

Ip = - {Zp + CIp- [G™'"] . (46) 

Substituting the expression for the port current into the aper¬ 
ture equations (l43l l yields an equation of the form of ( fT31 l for 
the aperture voltages 

ycav''j , ^ 2 /™'= . ( 47 ) 


Here the modified cavity admittance 


yrcav _ y^cai;_ 


R 


rad 


rad"] 


{Zp -f 


Up-Cap- [G' 


^adl ^/2 


(48) 

includes the effect of the current induced in the port antenna 
on the magnetic fields at the aperture. In the high loss limit 
{a > 1) or if the antenna load impedance is large this last 
term is small, r; amount of power 

coupled through the aperture is unaffected by the presence 
of the antenna. Once (E3 is inverted to And the aperture 
voltages, Va, dull can be used to And the current induced 
in the antenna, and the relation Vp = —Zpip can be used to 
And the port voltage. The power dissipated in the load is then 


Pp=Rp\Ip\^/2. 


We flrst consider the statistical properties of the fluctuating 
port voltage Vp in the high loss limit. As mentioned, in this 
case Z™"' ~ Z“’', and further ~ Z™^'- Thus, the 

aperture voltages Va determined by (l47l) will be the same as 
those determined in the radiation case of (fTTl) and the power 
coupled into the cavity will have the character displayed in Fig. 
|6] Statistical variations in the port current are then determined 
by the column vector ^ap in (l46l l. (In the high loss limit, 
variations in become small, and we can replace Z‘^°‘^ 

with for the port antenna.) 

The statistical properties of the elements of Uap the 
same as those of the off-diagonal elements of the general ma¬ 
trix ^ defined in (fl^ . In the limit of high loss the elements are 
complex, with independent real and imaginary parts, each of 
which are zero mean independent Gaussian random variables 
0, 0. The common variance for the real and imaginary parts 
is (27ra)~^. Since according to (l46T l the port current Ip is 
a linear superposition of Gaussian random variables, it too 
has independent real and imaginary parts that are zero mean 
Gaussian random variables. What remains then is to calculate 
its variance. Specifically, we And 



TDrad 1 

-o —Z4 • G™'' • Zi . 

\Zp + Z^‘^‘^\ Tra - 


(49) 


Similar arguments determine the statistics of the port volt¬ 
age. Since the port voltage and current are complex with 
independent Gaussian distributed real and imaginary parts, 
their magnitudes will be Rayleigh distributed with a mean 
determined by ( |49] |. Finally, we calculate the expected power 
coupled to the load (Pp) — Rp /2, 


(Pl) 


RpR^’^'^ 1 ^ 

\Zp + y™d|2 2-Ka ‘ 


(50) 


where Pt is the power transmitted through the aperture in the 
radiation case. We note that the ratio of the power coupled to 
the load and the power entering the cavity satisfies a relation 
identical to what we found previously for the case of localized 
ports HOl . 

The numerical computation of the exact expressions (|46]) . 
(l47l) . and (l48T l is now carried out. We now describe numerical 
solutions of the coupled port and aperture equations (I43II47I) . 
Elements of the matrix vector Up’ ^nd scalar f 

were generated using the same Monte Carlo algorithm used 
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Frequency [GHz] 


Fig. 10. Port power received by an antenna in the frequency range from 
1 GHz to 1.2 GHz. The cavity is highly irregular and overmoded, with loss 
factor of a = 1.0. The aperture is of dimensions L = 0.25 m xW = 0.02 
m, and the external plane-wave we assumed has amplitude = 1 mA/m, 

and direction of incidence 6 = n/4, -i/j = 0, and 0 = 0. Thin solid (black), 
dashed (blue), and dash-dotted (purple) lines indicate the port power for 3 
independent cavity realizations, while the thick solid (red) line indicates the 
ensemble average of the port power over 800 cavity realizations. Each cavity 
response is given by the superposition of 600 ergodic eigenmodes. 



Fig. 11. Port power received by an antenna in the frequency range from 
1 GHz to 1.2 GHz. The cavity is highly irregular and overmoded, with loss 
factor of a = 0.1. The aperture is of dimensions L = 0.25 m xW = 0.02 
m, and the external plane-wave we assumed has amplitude = 1 mA/m, 

and direction of incidence 6 = n/4, -i/j = 0, and 0 = 0. Thin solid (black), 
dashed (blue), and dash-dotted (purple) lines indicate the port power for 3 
independent cavity realizations, while the thick solid (red) line indicates the 
ensemble average of the port power over 800 cavity realizations. Each cavity 
response is given by the superposition of 600 ergodic eigenmodes. 


to produce Figs. (|7]l and Specifically, we generated 800 
realizations each with 600 modes. The cavity was excited 
through a narrow aperture of dimensions L = 0.25 m 
xW — 0.02 m, by a plane wave with amplitude |h*"“| = 1 
mA/m, direction of incidence 6* = 0, and ip — 7r/2, with 
polarization (j) = tt/2, and with frequency swept from 1 GHz 
to 1.2 GHz. The receiving port was an antenna with free-space 
impedance = 30 — j 20 H terminated by a Z/, = 50 H 
load. Figures [TO] and [TT] show the power coupled to the load 
for the case of loss factor a = 1.0 and a = 0.1 respectively. 
As in the cases of the power coupled into the cavity. Figs. 
|2| and |8| the power to the load for individual realizations 
shows variations with frequency which are more pronounced 
in the low-loss case than in the moderate loss case. Figures 
[Tol and [TT] also display the power to the load averaged over 
the 800 realizations. The amount of power reaching the load 
is generally in agreement with (fSOl) . From Figs. [7] and [8] we 
see that the average power entering the cavity at 1.1 GHz is 
1.6 X 10“® W and 0.6 x 10“® W in the a = 1.0 and a = 0.1 
cases respectively. Equation ( fSOl l then predicts for the power 
reaching the load 0.5 x 10“^ W and 2.1 x 10“^ W in these 
cases. Note that more power enters the cavity in the higher loss 
case, but more power reaches the load in the lower loss case. 
In the low loss limit we can develop a formula analogous to 
(i28T l. We assume the elements of the ^-matrix are dominated 
by a single cavity mode, and make the same approximations 
as the one leading to ( i28l l. The results for the power coupled 
to the load is found to be 


Pl{uj)=Pa — 
\ICo 


4:C(ACtP 

/C„ + * (ct + (Xa + Qfp)| 


(51) 


where Pa is the power coupled through the aperture at 
resonance, a a is the aperture loss factor defined in (i29l l. ap 
is an analogously defined loss factor for the port 


ap 




(52) 


and /Co = /C„ determines the shifted resonant frequency of 
mode-n, where 

/C" = K.n- 26 'aA-ap + Xl + R^‘^%p) /Rl ■ (53) 

In (fSTT i the statistically fluctuating quantities ap and a a 
determine the peak power at resonance reaching the antenna. 
In the low loss case considered here, the peak power delivered 
to the load can be a substantial fraction of the power passing 
through the aperture at the aperture resonance. Pa- 

V. Conclusion 

We have developed a statistical model of the coupling 
of electromagnetic power through an aperture into a wave 
chaotic cavity. The formulation combines the deterministic, 
and system specific elements of the aperture and a receiving 
antenna in the cavity, with a statistical model for the modes 
of the cavity. Particular attention has been focused on the case 
of a rectangular aperture, which has a set of well separated 
resonances. Power coupled through an aperture into free-space 
has peaks at frequencies corresponding to the modes of the 
aperture. When the aperture is backed by a cavity a new set 
of peaks appears at frequencies of modes of the cavity. These 
peaks can be as large as the free-space aperture resonance 
peaks. Simple formulas are derived that augment the detailed 
formulation, and that apply in the high-loss, low-loss, isolated 
resonance limit. These results are of interest for studying the 
coupling of an external radiation through apertures in complex 
cavities, such as RCs, and for practical scenarios of slotted 
enclosures populated by materials and electronic circuitry. 
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Appendix 


The admittance of a cavity excited through an aperture can 
be expressed by expanding the fields inside the cavity in a 
basis of electric and magnetic modes 

E = E (*) ’ (54) 

n 

and 

H = ^ (/^hr (x) + (x)) • (55) 


Here the electromagnetic modes satisfy the pair of equations, 
—iknC.^ = V X h®"*, and ifcnh®™ = V x e®'" with the 
tangential components of the electric field equal to zero on 
the cavity boundary including the aperture. The magnetostatic 
modes are irrotational h™® = —Vxn, where the potential 
satisfies the Helmholtz equation (V^ + k^) Xn = 0, with 
Neumann boundary conditions, |n- Vxnl^ = 0. The mag¬ 
netostatic modes are needed to represent magnetic fields that 
have nonvanishing normal components at the aperture BTll . It 
can be shown that all the magnetic field modes are orthogonal. 

The mode amplitudes are determined by projecting 
Maxwells equations onto the basis functions for each field 
type. The result of this action is an expression for the magnetic 
field amplitudes in the aperture that is equivalent to Q except 
that the radiation admittance matrix is replaced by a cavity 
admittance matrix 




■ j „,,em„,,em • „,,ms„,,ms 

tko I 


kl-kl V- 


ko 14"*® 


where 


= 

^sn 


/ dxj_es (xj.) • z X , 

Japerture 


(56) 


(57) 


is the projection of the magnetic filed of the cavity mode onto 
the aperture field profile and 


!/(■) 




(58) 


is a normalization factor for the eigenfunctions. 

Expression (l56l l is general and gives an exact expression 
for the admittance matrix of a lossless cavity in terms of the 
cavity modes. If we apply the random coupling hypothesis, 
we replace the exact eigenmodes with modes corresponding 
to random superpositions of plane waves. Specifically, near 
the plane z = 0 we write for the components of the electro¬ 
magnetic eigenmodes transverse to the z direction 

2 ^ 

h®™ = lim -^= ^ bj_L cos (kj • fiz) cos {6j + • xj_) , 

(59) 

where 9j are uniformly distributed in the interval [0,27r], 
|kjj = kn, with the direction of kj uniformly distributed over 
the half solid angle corresponding to kj-fi > 0, |bj| = 1, with 
bj uniformly distributed in angle in the plane perpendicular to 
kj. Except as mentioned, all random variables characterizing 
each plane wave are independent. A similar expression can be 
made for the scalar potential Xn generating the magnetostatic 
modes. With eigenfunctions expressed as a superposition of 


random plane waves each factor w^^i ) appearing in (l56l l be¬ 
comes a zero mean Guassian random variable. The correlation 
matrix between two such factors can then be evaluated by 
forming the product of two terms, averaging over the random 
variables parameterizing the eigenfunctions and taking the 
limit N ^ oo. We find for the electromagnetic modes the 
following expectation value 


\ ^ 


W 


(em) 



kl yk^ J 


(k X n) (k X n) 


(60) 


Here |k| = kn, represents the spherical solid angle of k, 
and V is the volume of the cavity. A similar analysis of the 
magnetostatic modes gives 


/ (ms) (n 

j Win W\ 


f d^k 


(k X n) (k X n) 

/ 27rE 

[Uv 

kl \ 


(61) 


The connection between the cavity case (l56]) and the radiation 
case (|4|i is now apparent. Specifically, we note that the factors 
are zero mean Gaussian random variables with a 
correlation matrix given by (l60l l. We can express the product 
in terms of uncorrelated zero mean, unit width 

s n 

Gaussian random variables by diagonalizing the correlation 
matrix. We again introduce matrix notation and represent the 
ss element of the product 


(em) (em) 

win W\ 

_ s n 

j/em 



1/2 


T r j^rad) ^/^ ) 

• WnWn • IS J \ , 

J ss 


(62) 


where is the radiation admittance matrix (|7]), and AA:„ = 
7r^/(fc^V) is the mean separation between resonant wave 
numbers for electromagnetic modes on a cavity of volume 
V. Substituting (l60l l into (l56l l we have 


-ycav _ 


2ikoAkn 


n (^0 - 


[g 


rad' 


1/2 


■WnWn 


j^radj 


+B^ 


(ko) , 
(63) 


where in the limit of a large cavity we have approximated the 
sum of the pairs of Gaussian random variables representing 
the magnetostatic contribution to the cavity admittance by 
their average values and using AA:„ = 27r^/ (Efc^) for 
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magnetostatic modes converted the sum to an integral 



i 


w 


{ms) 

sn 


W 


{ms) 


s n 


ko y™® 



(ms) (ms) 
Wsn 

s n 



k'^Akndflk 

~wr~ 


(k„_L X n) (k„j_ X n) 

L2 


iB'^ 


(ko) ■ 
(64) 


Equation (l6^ is the random coupling model prediction for 
the cavity admittance. The last steps are to replace the exact 
spectrum of eigenvalues, by a spectrum produced by 
random matrix theory and to insert a loss term. We introduce 
a reference wave frequency and associated wave number 
uJref = krefC, and we assume that the cavity is filled with 
a uniform dielectric with loss tangent ts = Under these 

assumptions for frequencies close to the reference frequency, 
the frequency dependent fraction appearing in (l6Tt can be 
expressed as 


where 


2fcoAfc 

Ui _ jL2 

«,g 


= [/Co - /C„ + ia\ 


-1 


^0 ^ ~ ^ref 

2fcoAfc Aw 


(65) 

( 66 ) 


measures the deviation in frequency from the reference fre¬ 
quency, and Aw = Akc, is the mean spacing in resonant 
frequencies. The resonant wave numbers are now represented 
as a set of dimensionless values 


i -2 _ 1-2 

^ref 

2koAk 


(67) 


which by their dehnition have mean spacing of unity. These 
are then taken to be the eigenvalues of a random matrix from 
the Gaussian Orthogonal Ensemble normalized to have mean 
spacing unity. Einally, the loss factor a is dehned to be 


_ 1 ^0 
““ 2^Ak 


( 68 ) 


We note that the zeros of the denominator are given by kg = 
kn — ia, which implies 5 (w) = —aAuj = —ojg/{2Q). Thus 
a = ujqI{ 2QAuj). This results in an expression analogous to 
those obtained in 0,0, m for the impedance matrix 


= i Im {X™''} -p ^ , (69) 


Thus, we have seen that in cases of ports described by planar 
apertures, we can express the model cavity admittance in terms 
of the corresponding radiation impedance or admittance and a 
universal statistical matrix 
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